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O ' Abstract 

(N 

j^' We investigate a minimal extension of the standard model in which the only new 

ingredient is the sterile neutrinos. We do not introduce extra Higgs multiplets or high di- 
mensional effective operators to induce mass terms for the active neutrinos, and the model 
is renormalizable in itself. We show for arbitrary numbers of generations and sterile neu- 
^ ■ trinos that the independent physical parameters in the leptonic sector are much less than 

previously anticipated. For instance, with three active and two sterile neutrinos, there 
are four mixing angles and three CP phases in addition to four non-vanishing neutrino 
^ ■ masses. We study phenomenological implications for tritium beta decay, neutrinoless 

O ■ double beta decay and neutrino oscillations. For the most natural see-saw parameters. 
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Q ■ we find that it is difficult to accommodate in the model the best-fit values of masses and 

r-| ■ mixing parameters from oscillation data no matter whether we include or not the null 

'p'' short-baseline experiments together with the LSND result. This implies that if the LSND 

^' result is confirmed by MiniBooNE, the see-saw parameter region of the model with two 

^ ■ sterile neutrinos could be largely excluded. 
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1 Introduction 

The experiments on neutrinos have offered the first piece of evidence that points to physics 
beyond the standard model (SM) of electroweak interactions [1]. The deficits in the solar 
electron neutrino flux [2] and in the ratio of the atmospheric muon to electron neutrino 
fluxes [3] can be best and most naturally interpreted in terms of neutrino oscillations [4]. 
These observations have been confirmed by experiments at accelerators and reactors [5, 6]. 
They imply that the three neutrinos are non-degenerate and interact non-diagonally with 
the charged leptons. Nevertheless, the controversial circumstance in the short baseline 
(SBL) experiments seems to indicate that the complete picture of the neutrino sector may 
be richer than with the three ordinary neutrinos: while the LSND experiment claimed 
to observe a statistically significant signal [7] that cannot be accommodated in the three 
neutrino scheme of solar and atmospheric experiments due to three very different mass 
gaps, the signal was not confirmed by other SBL experiments [8, 9]. This situation will 
hopefully be clarified in the near future by the MiniBooNE experiment [10]. 

The controversy at SBL experiments has stimulated a lot of theoretical attempts. 
Some of them appeal to certain drastically new physics like CPT violation [11], CPT 
violating quantum decoherence [12], or extra dimensions [13], to mention a few. A more 
conventional approach is to introduce additional neutrinos to accommodate the observed 
at least three mass gaps. Since the number of SM-like neutrinos is severely constrained 
by electroweak measurements to be three, these new neutrinos must be neutral with 
respect to the SM gauge groups; i.e., they must be sterile neutrinos. The concept of 
sterile neutrinos was introduced earlier in the other context [14]. The most economical 
scenario with one sterile neutrino has been extensively studied in the literature [15]. It is 
generally difficult, if not impossible, to explain all data in this scenario either because of 
the rejection of significant involvement of a sterile neutrino by the solar and atmospheric 
data or because of the tension existing between the positive SBL result at LSND and 
the negative one at all others. The next simplest would be to add two sterile neutrinos, 
the so-called (3+2) scenario. Sorel et al. have assessed carefully the compatibility of all 
SBL experiments, and found that the (3+2) scenario fits the SBL data significantly better 
than the one with a single sterile neutrino [16]. From the model building point of view, 
although it is simplest to add one sterile neutrino, there are viable models that contain 
two or three sterile neutrinos [17]. It has to be left for experiments to decide which is 
actually realized in Nature. 

In this paper, we consider a minimal extension of SM that could potentially accommo- 
date the neutrino data. It is minimal in the sense that only the neutrino sector is extended 
by adding some sterile neutrinos. We do not introduce additional Higgs particles or higher 
dimensional effective operators from some high energy scale. In particular, there are no 
bare masses for the ordinary active neutrinos. The model so extended preserves renor- 
malizability in itself. In an earlier work [18] we found that such a model has a strictly 
constrained leptonic sector which contains much less physical parameters than previously 



anticipated [19, 20]. With one sterile neutrino, for instance, there are only two mixing 
angles in addition to two non-vanishing neutrino masses, and the leptonic sector preserves 
CP symmetry automatically. Furthermore, the ratio of the two masses also appears in the 
mixing matrices, which makes the two otherwise independent mixing angles less effective 
in generating an experimentally favored mixing matrix than they would have been. The 
resulting model is even incapable of interpreting the solar and atmospheric data [18, 21]. 
In this work, we investigate the parametrization of the leptonic sector for any numbers 
of generations and sterile neutrinos. Since the number of physical parameters is much 
less than previously expected, it becomes numerically manipulatable to study its impli- 
cations on neutrino experiments even with two or three sterile neutrinos. For instance, 
with two sterile neutrinos, there are four mixing angles and three CP violating phases in 
addition to four non-vanishing neutrino masses, and the model could thus become viable 
to accommodate the neutrino data including the SBL experiments. 

The paper is organized as follows. In the next section, we study the leptonic sector of 
the minimally extended model, and count in particular the number of physical parameters 
contained in it for any numbers of generations and sterile neutrinos. To get prepared for 
phenomenological analyses, we work out analytically in section 3 the neutrino spectrum 
and leptonic mixing matrices for the most natural see-saw region in parameter space [22] . 
Then, in section 4 we discuss the phenomenological implications in the case of two sterile 
neutrinos. We do not attempt here a sophisticated statistical assessment; instead, with 
the fitting result by Sorel et al. for the general (3-1-2) scenario in mind, we consider 
whether it is possible to accommodate their result in the above mentioned parameter 
region of our minimal model. The answer turns out to be negative and even worse: even 
if we ignore the tension between the positive and negative SBL experiments, it is difficult 
to include the LSND result in the see-saw region of the model. We summarize our results 
in the last section and mention briefly further work worthy to do with the model. 

2 Mixing matrices for any numbers of active and 
sterile neutrinos 

In this section we first describe the leptonic sector of the n-generation SM extended by 
riQ sterile neutrinos. This is the minimal framework that can accommodate neutrino mass 
and mixing while preserving renormalizability of SM. By standardizing the neutrino mass 
matrices, we count independent physical parameters contained in the leptonic sector. 
Diagonalization of the leptonic mass matrices then yields the mixing matrices in the 
charged and neutral current interactions. For n > no, we show that the number of 
independent physical parameters can be further reduced due to the appearance of massless 
neutrinos. 



2.1 Setup of the model 

The only new fields compared to SM are the uq sterile neutrinos that we choose to be 
right-handed without loss of generality, srx, x = 1,- ■ ■ ,nQ. The model contains as usual 
the n generations of the lepton doublets, F^a = (^ia, fiaYi and of the charged lepton 
singlets, //?„, a = 1? ■ ' ■5'^- Here L,R refer to the left- and right-handed projections of 
the fields. 

Since the sterile neutrinos are neutral under SU{2)l x U{1)y by definition, they are 
allowed to have bare mass terms of Majorana type, 

I i 

^^SR = 2^^y^Rx^Ry + 2^*xy^^Rx (1) 

where ip^ = C'-f^ip* stands for the charge-conjugate field of ip with C = i'-f^'-f'^ satisfying 
C = — Cl" = —C^ = —C~^ and C'-f^^C = '-f^. We denote s^ = (si?)*" for brevity. The 
Uq X uq complex matrix M is symmetric due to anticommutativity of fermion fields, but 
is otherwise general. The presence of sterile neutrinos also introduces the mixing mass 
terms between active and sterile neutrinos through the Yukawa interactions, 

-^Y = yibFLaVfrn + y^xFLa^SRx + h.C. (2) 

where (p is the Higgs doublet field that develops a vacuum expectation value, ((/?) = 
(0, l)'^v/\^, and (^ = ia'^(f*. Note that no bare mass terms are allowed for the active 
neutrinos in this minimal extension of SM. The lepton mass terms are summarized by 



hmffn + nZDsR + h.c.J + - [s^rMsr + h.c.J (3) 



where rrif = y^v/\^ and D = y"'vl \pl are nxn and nxriQ complex matrices respectively. 



2.2 Standardization and diagonalization of mass matrices 

To facilitate counting the independent physical parameters contained in the leptonic mix- 
ing matrices after diagonalizing the lepton mass matrices, we first convert the neutrino 
mass matrices into a standard form. Since M is symmetric, it can be diagonalized by a uni- 
tary transformation sr — > YqSr to the real positive Y^MYq = Mdiag = diag(ri, ■ ■ ■ ,r„Q). 
The only other change in the total Lagrangian occurs in C^: D -^ DYq. Denote 
DYq = (oi, ■ ■ ■ , ttno) where a^ {x = 1, ■ ■ ■ , Uq) are general n-component complex vectors. 
Making a unitary transformation ul -^ X'^hl, we convert DYq into XDYo = D^ which 
has a zero triangle at its upper-left corner, and when n > no an additional (n — no) x no 
zero rectangle over the triangle. To keep the partnership of n/, and fi, we also transform 
/ -^ X'^ f so that the only other change is again in Cm which now becomes 






fiXnifX^fR + WDaSr + h.c] + - [sgMdiagS/j + h.c] (4) 



The matrix Da is found successively and in a way applicable to both cases n > uq 
and n < hq. First, we choose the unitary matrix Xi so that Xiai = (0, ■ ■ ■ , 0, |ai|)^ with 



ki| = y(i[(ii- Invariance of the inner product of vectors implies that {Xia^ 
determines all entries in the n-th row of Xi(ai, ■ ■ ■ , 



^^^ which 



l«i| 



). The first (n — 1) components of 
the column vector Xia^ are normalized to \l aJxPia^ with the projector Pi = 1 — ^^^. Next, 

) 



a^ai 



; '^no 



we choose the unitary matrix X2 that leaves all entries in the n-th row of Xi{ai, 
untouched and that rotates the first (n— 1) components oiXia2 to its (n— l)-th component. 
Invariance of the inner product then determines the (n — l)-th row of X2Xi{ai, ■ ■ ■ , a„g), 
and so on. The procedure continues by induction until the nmin-th column vector where 
ninin = min(n, no). Defining a^ = and the projectors 



which have the properties 

the end result, Z^a = X(ai, 02, • • 
vanishing entries. 
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• Xi, has the following non- 



(7) 



^n+l— a 



where the indices are restricted to a G [1, n] and a; G [n + 1 — a, no] when n < no, and to 
X G [1, no] and a E [n + 1 — x,n] when n > no. 

Now we diagonalize the lepton mass matrices. The charged part is done as usual by 
the bi-unitary transformations fi^ji = Xi^jii^ji with X'^{XmfX'^)Xji = mi being real 
positive. The matrix X\ then appears in the charged current interactions of leptons. To 
diagonalize the neutral part, we first rewrite the Lagrangian in terms of the fields (^^ j 

and their charge conjugates \c]- Then the neutrino mass terms become 
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where the (n + no) -dimensional, symmetric mass matrix in the new basis is 

/ 0, 



rrir. 



Dl 



Da 



(8) 



(9) 



with On being the zero matrix of n dimensions. Finally, the matrix m„ is diagonalized by 
the unitary transformation 

(i)-'-"^ (10) 

such that 

Y^ninY = m^ (11) 

is real and non-negative. The n x [n + no) submatrix, y, composed of the first n rows of 
Y appears in the charged and neutral current interactions of leptons. We shall denote the 
mass eigenstate fields of the neutral leptons, u, by the Latin indices, j,k = 1,2,- ■■ , u + uq, 
and those of the charged leptons, i, by the Greek indices, a, (3 = 1, ■ ■ ■ , n. Introducing the 
Majorana neutrino fields u = i^R + iyji, the changes in the total Lagrangian are summarized 
as follows. 
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where upon redefining X\ 



2.3 Counting of physical parameters 



(13) 



We note first of all that D/^ has different zero textures according to n > no or n < ng. 
This affects the number of physical parameters in the model. We start with the easier 
case of n < no. The matrix D/^ has n real entries lying at the rows and columns fulfilling 
a + a; = n + 1 and |n(2no — n — 1) complex ones to the right of the real entries. There 
are thus (n + no) real parameters and |n(2no — n — 1) complex ones in n2„, which are 
traded for (n + no) neutrino masses, |n(2no — n — 1) mixing angles and |n(2no — n — 1) 
CP violation phases in Y . Since for general original matrices M, D there are no further 
unitary transformations on sr^ til that leave the structure of n2„ invariant, all those 
angles and phases must appear in the submatrix y of F, and thus in V^ and V'^ . The 
additional parameters introduced to V'^ on diagonalizing the charged lepton masses are 
counted as usual: \n{n — 1) angles and ^n{n — 1) phases after absorbing the n phases 
by the charged lepton fields. Again, for general original matrices, they are not expected 
to combine with or cancel those in y. In summary, there are n massive charged leptons, 
(n+no) massive neutrinos, n(no — 1) mixing angles and n(no — 1) CP phases in the charged 
current matrix V^ . Out of those, only in(2no — n — 1) angles and |n(2no — n — 1) phases 
appear in the neutral current matrix V^ of the neutrinos. 

When n > no, in addition to the indicated zero triangle, the first (n — no) rows of 
D/x also vanish, which correspond to (n — no) massless neutrinos after diagonalization. 



Arbitrary unitary transformations, yo, amongst the massless modes are allowed without 
changing the matrix ?72„. Thus, y has the structure 
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(14) 



While both y and yo appear in V , only y appears in V . In this case, there are uq real 
and 1^0(^0 ~ 1) complex parameters in Z^a; together with Mdiag, they are traded for 2no 
neutrino masses, |no(no — 1) mixing angles and |no(^o — 1) CP phases in y and thus V^ . 
Now we count the physical parameters contained in V'^. An n dimensional unitary 
matrix V may be parameterized as a product in any arbitrarily specified order of the 
n diagonal phase matrices, ea{ua) (« = l,---,n), and the ^n{n — 1) complex rotation 
matrices in the (q;,/3) plane, ujap{(^ai3,<^ai3) {n > j3 > a > 1) [19]. Here ea{z) is obtained 
by replacing the a-th entry in 1„ by the phase z, and 



0Oa/3{0^(3,Va/3) = e„(e'^«^)i?,^(^„^)e,(e-''^«'3) 



(15) 



where RapiOap) is the usual real rotation matrix through angle 6aj3 in the (a,/?) plane. 
We choose the order of products in such a way that it fits our purpose here: 



V = V^V^Vo 



(16) 



where Vq is the general unitary matrix in the subspace spanned by the first (n — no) 
axes, V2 the unitary matrix in the subspace spanned by the last n^ axes, and Vi is the one 
mixing the two subspaces. In V^ ^ Vq can be cancelled by y^ in y*, which remains free until 
now. The no diagonal phases in V2 can be arranged to its very left to get absorbed by the 
charged lepton fields. We are thus left in V'^ with the |no(no — 1) complex rotations in V2 
and (n — no)no ones in V^i, in addition to the |no(no — 1) complex parameters contained 
in y. 

However, not all parameters in Vi are physical. To see this, we write 



Vi = YluJaz{6az,Va 



(17) 



with a E [l,n — no], z E [n — no + l,n]. The order of products is specified as follows: 
factors with the same z are grouped together with z increasing from left to right; within 
each group, a also increases from left to right. Denoting for brevity ujazi^az, Vaz) = ^az) 
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(19) 



The left phase matrix in the first hne can be pushed through V2 to be absorbed by the 
charged lepton fields, while the right phase matrix in the last line passes through y in 



which ?/o 



1 



n— no 



to get absorbed by the massless neutrino fields. This leaves with us 



the (n — no)(^o ~ 1) phase differences and all of the real rotations in Vi that cannot be 
removed. To summarize, there are no(n — 1) physical mixing angles and n(no — 1) physical 
CP phases in V'^. 

There is another way to count the physical parameters in V'^. A general n x {n + Uq) 
complex matrix V'" has 2n{n + no) real parameters. First, V'^V'^'' = In offers n^ real 
constraints. Second, V*^ may be multiplied from the left by an arbitrary diagonal phase 
matrix by redefining the charged lepton fields without affecting physics. This amounts to 
another n real parameters. Finally, since there is no bare mass for the active neutrinos, 
we have y*m^y'^ = 0„ which implies V'^m^V'"'^ = 0„ [23]. The constraint is symmetric, 
and thus imposes 2n real conditions from the diagonal and n{n — 1) from the off-diagonal. 
Thus, the number of real physical parameters in V'^ is 2n(no — 1). However, this counting 
is not as advantageous as the above one. It does not distinguish between mixing angles 
and CP phases, or tell how many of them enter in V^ . More importantly, it does not 
take into account the further reduction of physical parameters due to the appearance of 
massless neutrinos when n > tiq. 

The counting of physical parameters is summarized in Table 1 where the result of 
Ref. [19] is also shown for comparison. The difference arises from the fact that the 
zero bare mass for active neutrinos has been completely exploited here to remove all 
unphysical parameters while it was only partially applied in Ref. [19] to delete unitary 
transformations within the massless neutrinos for the case n > Uq. 



3 Approximate results for two sterile neutrinos 

The number of active neutrinos has been constrained by experiments to be n = 3. To 
accommodate neutrino mass and mixing, one would introduce sterile neutrinos as few as 
possible. For no = 1, we have two massless neutrinos and two massive ones; and there are 
two mixing angles in V'". This has been shown in Refs. [18, 21] to be even impossible to 
explain the solar and atmospheric experiments that call for two mass squared differences 
(Am^j = rrij — mf) and two mixing angles, let alone the LSND results hinting at a third 
Am^j. The next simplest is the case with two sterile neutrinos which is studied in this and 



Table 1: The numbers of independent physical parameters are shown for the mixing matrices 
V^ and V^ . Note that all parameters in V'^ are aheady included in V'^. Also shown are the 
resuhsofRef. [19] for V^ 
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the next section. Now we have one massless neutrino (v\) and four massive ones (z>'2,...,5) 
Separating out the massless mode from the matrix Z^Aj we have 
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(20) 



where generally (ii 2 > 0, ri 2 > and z is complex. These parameters are traded for the 
four neutrino masses n22,...,5 and one mixing angle plus one CP phase in V^ . Both the 
angle and the phase enter in V^ which includes additional three angles and two phases 
coming from the matrices Vi, V2. We thus should have enough degrees of freedom to 
accommodate oscillation experiments that require at least three Am^j and four mixing 
angles. 

Although it is possible to diagonalize algebraically the mass matrix n2„, it is not in- 
structive for our phenomenological analysis. Experimentally, we need three well-separated 
Am^j to explain the positive oscillation results, Ami, ^ Am^TM ^ ^^lsnd- Since the 
massless neutrino vi is purely active and the solar and atmospheric experiments are more 
in favor of active neutrino mixing than the involvement of sterile neutrinos, the other two 
mainly active neutrinos z/2,3 have to lie close to vi while the two mainly sterile neutrinos 
z/45 must be well above it; i.e., m^^^ ^ ni2,3- Furthermore, since the sterile mass terms 
are not constrained by the low energy symmetries of SM, it is natural that they may 
be linked to some new physics at a higher scale. Thus, we shall assume the hierarchy 
^1,2 ^ di,2i \z\ in our phenomenological analysis. As a special case, we shall also consider 
the possibility of 2; = in which analytic results can be easily worked out for masses and 
mixing matrices without assuming the hierarchy ri 2 ^ (ii,2- 
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3.1 See-saw case: ri 2 ^ di^2, \z\ 

Excluding the massless neutrino, the mass matrix to be diagonahzed via Y'^rrinY 
diag(m2, ■ ■ -,^5) is 

/ d2 \ 

di z 
d\ ri 
\d2 z r2 J 

The eigenvalues exact to the second order in the expansion of (ii,2,-2 over ri^2 are found 
to be 



(22) 



^^2 3 = 


1 


dl 1 {z + d2p*^f 




1 + Ipt ' 


n r2 




, dl 




77^4 = 


= ri + ^ 





1715 



££, 



r2 + 



s + kp 



^2 



where ra2 < m^, 17121713 = -^-^ and 
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When z is real, there is an additional useful relation, 

d\ d\ + z^ 

7712+1713 = \ 

ri r2 

The first two rows of Y give the submatrix y in y. Including the massless degree of 
freedom corresponding to yo = 1 after removing unphysical parameters, we have 

/ 1 \ 

y*=\ V-P- V+P+ J (25) 

\Q v_ v+ ^ ^ / 

where the second and third columns are exact to the first order and the last two to the 
second order, and 
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where 2q;=p are the phases of the complex numbers. 
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Parameterizing the 3x3 unitary matrix V = (Vap) with q;,/5 = e,p,T, we have 
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Since y^y = I3 is exact at the first order in the expansion, so is V'^V'^'' = I3. 
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3.2 Case: z = 

The neutrino mass matrix can be easily diagonalized in this case. The masses are 
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There is no free angle or phase in the diagonalizing matrix Y which is completely fixed 
by the masses, so that 
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"12 +"1-5 V m2+m5 

Note that the i factors are introduced to make m2,3 real positive and do not by themselves 
signal CP violation [24]. Then 



Vee iVe^S2 iVerSl VerC\ Ve^C2 

V^ = I Vf,e iV^^S2 iV^rSl V^rCl V^^C2 

Vre iVrf,S2 iVrrSl KrCl Vrf^C2 



(31) 



We have \/^\/^'l" = I3 exactly. 



4 Implications on neutrino oscillations 

We consider now the phenomenological implications of the previous section. Before we 
move to neutrino oscillations, we discuss briefly the results on the tritium decay [25] and 
the neutrinoless double (3 decay [26]. Both decays are currently available experiments 
sensitive to the absolute neutrino mass. The first one measures the following effective 
mass through the distorted decay spectrum. 



m 



l = T.^'\VS\'={Vy*mly^V% 



(32) 



Using Y'^rUnY = m^, y*ra1y^ is just the upper-left 3x3 submatrix oininral^, i.e., D^D^^. 
Thus the effective mass is only sensitive to the Dirac mass terms. 



m„ 



dl\Vef,\^ + {dl + \z\^)\Ver\^ + d2{zV*Ver + Z*Ve^V, 



(33) 



which is of order the product of a light neutrino mass and a heavy one if Vee is not very 
close to unity. 

The neutrinoless double j3 decay violates the lepton number conservation and can occur 
only when the neutrino is of Majorana character. At the leading order in the expansion 



11 



of neutrino mass over the characteristic momentum transfer in the nuclear transition, the 
decay is proportional to the effective mass 

mee = lE^.(K5)'l (34) 

Note that there is interference amongst different terms. In the model considered here 
where the only extension is to add sterile neutrinos to SM, the above mass vanishes. 
Actually we can show that the effective mass for neutrinoless double like-sign charged 
leptons (^a^J) decays 

^a/3 = lE^^KS-^ftl (35) 

vanishes generally for any number of generations and any number of sterile neutrinos. 
Using V'^ = Vy*, we have 

n+no 

Y. ^J^aj^S = H y^^Vf^s Yl ^jy^jVlj = Y ya^yps{y*m^y^)^s (36) 

j=l 7,<5 j 7,(5 

where {y*iTLi,y'')^s constitutes the upper-left n x n submatrix of Y*mi,Y'' = rUn, which 
vanishes however due to the absence of Majorana-type mass terms for active neutrinos in 
the current model. While this does not necessarily mean that such decays are forbidden, 
because they can be induced at the next orders of the expansion and weak interactions, 
it does imply they are strongly suppressed. It could also be the case that one or more 
neutrinos are much heavier than the momentum transfer so that the expansion does not 
apply to them and the sum over j is not complete. Even in this case, the decays are 
still suppressed since the vanishing of the complete sum and very large masses of those 
neutrinos imply very small couplings between them and the charged leptons. 

In the subsections to follow, we will consider the potential of accommodating neutrino 
oscillation experiments in our simple model. If the LSND signal does exist at all, it must 
point to a much larger mass squared splitting than those in the solar and atmospheric 
experiments. Since the lightest active neutrino is massless in our model, it implies that 
the masses of the neutrinos responsible for the LSND signal are much larger than those 
of mainly active neutrinos involved in the solar and atmospheric experiments. It is thus 
adequate to work in the see-saw hmit (ri,2 ^ '^1,2, l-^l) to simplify our numerical analysis. 
We will also consider a special case with z = in which analytic results have been readily 
worked out in the last section. 

4.1 See-saw case: di,d2, \z\ <C ri,r2 

4.1.1 Including LSND and negative SBL results 

A difficulty to accommodate all oscillation results by adding sterile neutrinos is the tension 
existing between the negative and positive SBL experiments themselves. The statistical 
compatibility of the two sides has been analyzed in Ref. [16] amongst others [27]. It 
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was found in Ref. [16] that the SBL results can be significantly better reconciled with 
two sterile neutrinos than with a single one. Assuming the SBL compatibility and CP 
conservation, they found the following best-fit point with two sterile neutrinos, in our 
notation. 



rriA 



0.92 eV^ 
= 0.121, 



m? 



22 eV' 
0.204, 



IK^I = 0.121, |V;C| = 0.204, \VM\ = 0.036, \V^,\ = 0.224 
Restricting the neutrino masses to sub-eV, the best-fit point moves to 



(2) 



ml = 



0.46 eV^ 
= 0.090, 






0.89 eV^ 
0.226 



vM\ 



0.125, I ye I 



0.160 



(37) 



(3^ 



The above results were obtained for a general model of three active and two sterile neu- 
trinos, in particular without parametric constraints on the mixing matrix V'^. In our 
minimal model, the independent physical parameters in V'^ are much less than the gen- 
eral case as detailed in the last section; and in addition, V'" involves the neutrino mass 
ratios. It is thus interesting to ask whether the above best-fit points can be accommodated 
with the minimal number of physical parameters in the model. 

With two sterile neutrinos, there are four mixing angles in V'" = Vy*, three of them 
from V and one from y*. Assuming CP conservation, z and V become real so that no 
phase can appear. Although the second and third columns of y* are pure imaginary, their 
i factors have no effect on oscillations. Since the above best-fit points fall in the see-saw 
limit, it is good enough to identify m^ ^ ri, m^ ^ r2 without loss of generality. For the 
solar and atmospheric mass gaps that enter our analysis through the matrix y*, we use 



Ami 



8 X 10"^ eV 



^"^ATM 



2.4 X 10"^ eV 



(39) 



A slight deviation from those numbers will not change our qualitative conclusion. Since 
ms > 7712 > mil = 0, we must have 7n\ k, Am^rpj^ while there are two possible ways to 
arrange for the solar gap: either (a) 7n\ ^ Am^ or (b) ^711^2 ~ AttIq. Combining with 
the two best-fit points shown in eqns. (37,38), we have four schemes to consider, named 
below as (la), (lb), etc. 

We find that schemes (la,lb,2b) can be excluded without using any constraints from 
the solar and atmospheric mixing angles. Take scheme (la) for instance. The best-fit 
values for the matrix entries in eqn. (37) translate into 



v;,|J = 0.121, 

^ + V —I 



w. 



^ - 0.204, 
0.036, \V..^ + V.r- 



\V |2i 



0.224 



(40) 



The first two give ^ > VO.1212 + 0.204^ ^ o.237. Using eqn. (24) we have the bound. 



< 



/m2 — a/ttis)^, which implies then — < 



I v^m?- \/ma"| 



5.85 X 10 ^. Now we multiply 



the last one in eqn. (40) by ^, use the second and 
at the following result 



did2 
rir2 



m2m:i 
m4m5 



0.987 X 10"V„„ ± 0.204^ 

r2 



0.224— 

n 



0.987 X 10"^ to arrive 



(41) 
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where the sign in the second term refers to V^r^- Even if the second term takes its largest 
magnitude of 1.20 x 10~^ and adds to the first one, it still requires |V^^| ^ 1 to reach 
the lowest value of the right-hand side of 5.31 x 10~^. This scheme is thus excluded, and 
similarly with the schemes (lb, 2b). 

For the scheme (2a), we have to appeal to a constraint from the solar mixing angle. 
The best-fit matrix entries correspond to 



|Kr|J = 0.090, |V;.|J = 0.226, 

|K4 + Verf-J = 0.125, |V;4 + V,rf;\ = 0.160 



(42) 



The third one gives 



0.0262Ku ± 0.090^ 

^2 



0.125— (43) 



n 



where ^ > 0.243, ^ < 0.131 and we have used '^^ = 2.62 x lO^^. Since the best-fit 

C _ T/ ~ . /l _ oivi2 f, „^^ rio,ro IT/ I ^ o _ oivi /3 ~ 



solar angle 9q is about 32°, using V^ = Vee ~ yl — sin 9q we have \Ve^\ < s© = sin^© 
0.53. Then the left-hand side in the above equation cannot exceed (1.39 + 1.18) x 10^^, 
significantly deviating from the lowest value of the right-hand side, 3.04 x 10~^. The 
difficulty can be better viewed the other way around. Assuming V^^ + V^^ = s| and using 

^ = !^^^±^ _ f^y !!^ < 1.896 ^ 10- (44) 

ri m^ \rij mg 



the LHS of the third equality in eqn. (42) cannot exceed the value, JV^^ + V^^J ^^^ < 

O.138s0. Then, the equality requires that Sq > ^^ = 0.908, or 9q > 65°, which deviates 
far away from the best-fit solar angle. 

4.1.2 Excluding negative SBL results 

Assuming the positive and negative SBL experimental results are compatible, we have 
seen in the above that the best-fit points obtained in Ref. [16] cannot be accommodated 
in the see-saw parameter region of the minimally extended model which has much less 
physical parameters than a general nonrenormalizble model. The deviations are so large 
that we suspect it may be difficult to find a parameter region to encompass all SBL results 
that is consistent with other oscillation experiments. Although this issue should better 
be assessed by an extensive statistical analysis, it goes beyond the scope of this work. 
Rather, we would ask the question: since it is generally difficult to relax the tension in 
SBL experiments, what will happen if we take seriously the LSND signal but ignore the 
negative SBL results? 

As the mass gap required by LSND is much higher than those by other experiments, 
a simple but plausible way to answer the question is to take the best-fit values for the 
latter and consider whether there is any room for the LSND signal. Since only LSND 
hints at a third mass gap besides the known solar and atmospheric ones, it is natural to 
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simplify our analysis further by arranging either ri <^ r2 or ri ?a r2. In the first case, 
the neutrino u^ essentially decouples whose only effect is to provide a small mass to z/2- 
The decoupling is verified by the small entries in the 5th column of V^. In the second 
case, z/4 5 effectively form a Dirac neutrino when m^ = m^ holds exactly. Even if they 
have a gap that happens to be a solar or an atmospheric one, their contributions to other 
experiments can still be ignored because of the smallness of the relevant mixing entries. 
We assume CP conservation so that V and z are real. For definiteness, we assume 
2; > so that p_|_ > 0,p_ = — p+^ < 0. The conclusion does not change for z < 0. For 
oscillations, we can also ignore i factors in the 2nd and 3rd columns of V'". Then all 
matrix entries are effectively real. The oscillation amplitude responsible for the LSND 
signal appropriate for both cases of ri <^ r2 and ri ^ r2 is 






(I + fj) KrV + I^^V;^ + ^{Ve^V^r + VerV^^] 



(45) 



which involves the upper-right 2x2 sub matrix in V^. To determine it, we apply the best-fit 
point {9q = 32°, ^atm = 45°, ^chooz = 0) to the relevant 2x2 submatrix in V*^: 

T/C _ V'erP+-VeM _ „ T/C _ Ve^p++VeT _ r. 

which solves the V^ entries in terms of 6q and p^. The amplitude becomes 



^LSND — 2S0 



'^ ^ f!^ P+(1 + C0P+) , dl Cq - p+ d2Z p\ - 2cqP+ - 1 



When ri ^ r2, we have approximately 



(47) 



1712 ~ — , ms ^ — , m4 ^ ri, p+ ^ —5 < 1 (4j 

r2 ri df r2 



The amplitude is dominated by the first term 



Alsnd - 244p+ - 2444 - 24^4 (49) 



which is very small because 2; <^ ri(~ 1 eV ) <^ r2 and m2 = JAttiq < 0.01 eV . 
For an order of magnitude estimate, we note that z should be about the same order 
of di2- Then, the order of Alsnd should be roughly from 2s% . 2 ® ~ 4 x 10^^ to 

2g2 /^P^l^ < 2 X 10-^ which is well below the required LSND level of ~ 3 x lO^^. 



For 1714 ^ ri ^ r2 ~ ^ArriLSND) ^^ have the relations 

p-' - p+ ^ ^'~fv ^ . ^2^4 ^ ^^^ ' ^ ^7 : :i "^^^"" (50) 
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_i rf? - rf^ + z'^ {dj + z^)p\ + dl- 2d2zp+ 

,^ -p+ ^ d^z ' ^^^^ ^ TT^ 



which, after some algebra, simphfy considerably the amplitude to 

in /K in 

which is again far below the desired level. 

In summary, even if we ignore the negative SBL results to avoid the potential in- 
compatibility in SBL experiments, we cannot naturally accommodate simultaneously the 
results of solar, atmospheric, reactor and accelerator experiments including LSND in this 
model with see-saw parameters di,d2,\z\ <^ ri, r2. 

4.2 Case: z = 

As we argued in the above, the see-saw hierarchy ri^2 ^ (^1,2, l-^l is natural both theoret- 
ically and phenomenologically. The conclusion reached in the last subsection is thus 
quite general. Nevertheless, as a good illustration we consider below a special case 

where all relevant quantities can be calculated without working perturbatively. Since 
5 

]^m^ = det(m|^m„) = (^1^2)"^, neither of rfi 2 can vanish. We therefore study the case of 

i=2 

2; = in a way parallel to the last subsection. 

4.2.1 Including LSND and negative SBL results 

Although rus < 1714, 1712 < m^ but otherwise arbitrary, theoretical and phenomenological 
considerations require that ri^2 ^ '^i,2- Thus we generally have 7724^5 ^ ^^2,3- As the 
problem is symmetric under the simultaneous interchanges 1^2 ^^ 1^3 and u^ ^^ u^, we can 
choose without loss of generality m^ > m^. Since the relative ordering of 7722,3 is still free, 
we have four schemes to arrange for the solar and atmospheric mass gaps: 

(51) 7T7^ = Att^I, ml = Aml^^; 

(52) ml^ = Ami, ^l = ^"^atm! ^^2) 

(53) ml = Ami, ^2 = ^"^atm! 

(54) 777^3 = A777I,, 777^ = A777iTM 

The parameters Cj, Sj are collected in Table 2 that are computed for the best-fit values 
shown in eqn. (39) and the best-fit values 7774, 7775 of Ref. [16] shown in eqs. (37,38). 
A slight shift of the best-fit values does not change much those parameters because of 
the double square root dependence, and in particular our qualitative conclusion below is 
rather stable. 

It is clear from the table that it is far from possible to get close to the large best-fit 
values for |V^|, |V^^| in eqn. (37). For the sub-eV fit, schemes (S2,S4) seem to have 
a chance, so let us take a closer look. We find |Ver| ~ 0.35, \Ve^\ ~ 0.56, \Vfj,r\ ~ 
0.87, \Vf,f,\ ^ 0.73. The unitarity of V is violated by \V^r\^ + |V;,^p ^ 1.3 > 1. More 
seriously, large \Vfj,r\ and \V^fj,\ imply a very small \Vfj,e\ = |V^| which cannot be tolerated 
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Table 2: The parameters Cj, Sj are computed for the best-fit values ArriQ = 8 x 10 
Am^rpj^ = 2.4 X 10~^ eV^ and the best-fit values ml, ml shown in eqs. (37, 38) 
sponding to the left and right part of the table. 



-5 eV^, 
corre- 



schemes 


Cl 


Sl 


C2 


S2 


Cl 


Sl 


C2 


S2 


(SI) 


0.22 


0.98 


0.04 


1.00 


0.26 


0.97 


0.10 


0.99 


(S2,S4) 


0.22 


0.98 


0.10 


0.99 


0.26 


0.97 


0.22 


0.98 


(S3) 


0.10 


0.99 


0.10 


0.99 


0.10 


0.99 


0.22 


0.98 



by the data of Am^rpj^ oscillations; the circumstance with |VeT-| and \Ve^\ is conflicting: 
while their intermediate size implies a too large \Vee\ that is excluded by Chooz, they are 
not large enough to explain Am^ oscillation data. We can thus safely conclude that the 
best-fit points of Ref. [16] cannot be realized in the current model. 



4.2.2 Excluding negative SBL results 

Now we ignore the negative SBL experiments and consider the consistency to accom- 
modate all other experiments in the current model. It is sufficient to restrict ourselves 
to those that can be reasonably well described by vacuum oscillations, i.e., KamLAND, 
K2K, Chooz, and LSND that cover all three mass gaps. We take as our input the mass 
gaps shown in eqn. (39) and AmLg^D ~ 1-2 eV^, and study the possibility to obtain 
desired oscillation amplitudes for the mentioned experiments. Again we consider the nat- 
ural case 7714^5 ^ ^2,3 so that 2/4^5 are mostly sterile and 1^2,3 mostly active, and simplify 



mt 



the matter further by working with the 'only three gaps', i.e., either m 

or ms ^ m,4 = Attilsno- The masses m,2,3 are classified into the four schemes shown in 

eqn. (52). 

For m,4 = m^, schemes (SI) and (S3) are equivalent because they are related by a 
renumbering of indices; similarly with (S2) and (S4). Thus we only need to consider the 
schemes (S1,S2) below. The masses yield 



(SI): 
(S2): 



Cl 
c = 



^ 0.21, si ^ 0.98, C2 ^ 0.09, 

Cl ^ C2 ^ 0.21, S = Si ^ S2 



S2 ^ 1.0 

^0.98 



(53) 



The oscillation amplitudes that are multiplied to the oscillating factors are, in scheme 
(SI): 



LSND 

Chooz 

K2K 

KamLAND 






e3l 






+ 



e2l J 

^^l2\ ) 



AW V* r^ 
A\y pfi _ IT/ 



V V* r^P 

2^„2 



'1 



MV, 



fJ.T\ 



'i-\v, 



/IT\ 



4\V r\V 

^\ '^ ee\ K e/i| 



(54) 
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and in scheme (S2): 



LSND 

Chooz 

K2K 

KamLAND 



AHI/Cp I |T/C|2\|T/C|2 ^ A\V 12(1 _\y \2\„2 



(55) 



12^4 



' fJ.T\ 



2 ^ 



Consider scheme (SI) first. Positive K2K and negative Chooz results imply |V^ 
|V^T-P ~ I and tiny iVg^P- The order 10"^ exclusion level at Chooz restricts the cf term 
for LSND to be less than ~ 10~^, thus subdominant to the C2 term, which however 
cannot exceed c| ~ 6.6 x 10"^, well below the desired LSND signal. Scheme (S2) is not 
better. K2K and Chooz imply |V^eP ~ | and tiny \Vee\^- But in that case, the LSND 
result should be more negative than the Chooz because the ratio of their amplitudes is, 
g^^^ Ri ^|V^eP ~ 10~^, in sharp conflict with the observations. 

When ms ^ 7724 ^ ^^2,3, we have C2 ~ 0, S2 ~ 1 and v^ decouples from the os- 
cillation. Since 7722 ~ — , it is more appropriate to assume m^ = Am|, than any other 
arrangements. Thus we only have to consider the scheme (SI). This is just a special case 
of the above discussions with the difference that the LSND signal becomes more difficult 
to accommodate since it arises from the first term that was ignored in the above. 

5 Conclusion 

If the LSND experiment is confirmed by MiniBooNE, the scenario of three ordinary active 
neutrinos will be insufficient to explain all neutrino oscillation data. A natural approach 
to the problem is to add sterile neutrinos to allow for more independent mass squared 
differences that could potentially accommodate three well-separated mass gaps. We have 
investigated systematically the simplest type of such models where the only extension 
to SM is the addition of sterile neutrinos. In particular, we do not introduce extra 
Higgs multiplets or higher dimensional effective operators to induce masses for the active 
neutrinos, and the extended model keeps renormalizable as SM. We found that for any 
numbers of generations and sterile neutrinos its leptonic sector contains much less physical 
parameters than previously expected, upon exploiting completely the texture zero in the 
neutrino mass matrix. The model thus becomes quite viable for phenomenological analysis 
even if it contains two or three sterile neutrinos. We demonstrated that the mixing matrix 
in the leptonic charged current interactions has a factorized form and that the factor 
containing the neutrino mass ratios makes the neutral current interactions of neutrinos 
non-diagonal as well. 

We have studied the phenomenological feasibility to accommodate all oscillation results 
in the extended model with two sterile neutrinos. We have restricted ourselves in this 
work to the see-saw region of the parameter space which is most natural considering 
the experimentally found well-separated mass gaps and for which all relevant results can 
be worked out analytically. We have used as our reference points the best-fit values of 
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Ref. [16] based on the analysis of the complete SBL data. Unfortunately, the answer 
is negative, and even more: even if we take the LSND result seriously but ignore other 
null SBL results, it is difficult to accommodate the best-fit values for the solar (plus 
KamLAND), atmospheric (plus K2K), Chooz and LSND results in the see-saw parameter 
region. A slight shift of the best-fit values does not alter our qualitative conclusion. We 
attribute this to a new feature of the minimally extended model exposed here: the leptonic 
charged current mixing matrix is strongly modified by mass ratios of neutrinos. Adding 
more free mixing angles does not necessarily improve the simultaneous fitting of masses 
and mixing angles. If the LSND signal is confirmed by MiniBooNE, the see-saw region of 
the model with two sterile neutrinos could largely be excluded. 

Finally we discuss briefiy how the current work could be extended. Although the see- 
saw parameters are most natural from theoretical point of view, it is highly desirable to 
make a complete scanning of the parameter space. Since the mainly sterile neutrinos are 
generally light or at least not very heavy in the model, fitting data in regions other than 
the see-saw one may not cause too serious fine-tuning of parameters. Suppose the LSND 
result will be confirmed by MiniBooNE, a negative result of the scanning would rule out 
the simplest extension of SM in its leptonic sector and call for something really new to 
SM. 
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